“-A134  276  A  GENERALIZED  COUPLE  THEORV  OF  HVGROTHERMflL-ELASTICfTV:  1/i 

TRANSIENT  EFFECTS.  .  (U)  LEHIGH  UNIV  BETHLEHEM  PA  INST  OF 
FRACTURE  AND  SOLID  MECHANICS.  .  G  C  SIH  ET  AL.  SEP  83 
UNCLASSIFIED  AMMRC-TR-82-56  DAAG46-82-K-0006  F/G  11/4  NL 


Q> 

Lad  .  I 


AMMRC  TR  83-56 

A  GENERALIZED  COUPLE  THEORY  OF  HYGROTHERMAL ELASTICITY: 
TRANSIENT  EFFECTS  IN  COMPOSITE  LAMINATE  WITH  CIRCULAR 
CAVITY 

September  1983 

G.  C.  SIH  and  M.  T.  SHIH 

Institute  of  Fracture  and  Solid  Mechanics 

Lehigh  University 

Bethlehem,  Pennsylvania  18015 

FINAL  REPORT  Contract  No.  DAAG46-82-K  0006 


Approved  for  public  release;  distribution  unlimited 


Prepared  for 

ARMY  MATERIALS  AND  MECHANICS  RESEARCH  CENTER 
Watertown,  Massachusetts  02172 


DTIC 

ELECT E 

NOV  3  1°83 


r 


D 


B 


83  it  03  095 


The  findings  in  this  report  are  not  to  be  construed  m  an  others* 
Department  of  the  Army  position,  unless  so  designated  by  othei 
authorised  documents 

Mer.uon  ohany  trade  names  01  manufacturers  in  this  report 
shall  not  be  construed  as  advertising  noi  as  an  official 
indorsement  or  approval  of  such  products  or  companies  by 
the  Uenrt  States  Government 


DISPOSITION  INSTRUCTIONS 


Destroy  thu  reunrt  wNn  it  is  no  longer  needed 
Do  not  return  it  to  the  originator 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  face  (WTiwi  Dtim  Emt'td) 


|  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

AMMRC  TR  83-56 

wmm 

\  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (* ml  Submit) 

A  GENERALIZED  COUPLE  THEORY  OF  HYGROTHERMAL- 
ELASTICITY :  TRANSIENT  EFFECTS  IN  COMPOSITE 

s  type  of  report  a  period  covereo 

Final  Report  - 
1  Nov  81  to  1  Jan  83 

LAMINATE  WITH  CIRCULAR  CAVITY 

6  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHOR!  •  > 

G.  C.  Sih  and  M.  T.  Shih 

t  CONTRACT  or  GRANT  NUMBER!*; 

DAAG46-82-K-0006 

9  PERFORMING  ORGANIZATION  name  AND  aooress 

Institute  of  Fracture  and  Solid  Mechanics 

Lehigh  University 

Bethlehem,  Pennsylvania  18015 

10  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  «  WORK  UNIT  NUMBERS 

D/A  Project:  8H363304D215 

AMCMS  Code:  603000.2156 

II.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

Army  Materials  and  Mechanics  Research  Center 

12  REPORT  OATE 

September  1983 

1  ATTN:  DRXMR-K 

I  Watertown,  Massachusetts  02172 

•)  NUMBER  OF  PAGES 

51 

1  M  MONITORING  AGENCY  NAME  A  AOORESSfif  different  from  Controlling  Ofhre) 

IS.  SECURITY  CLASS  (ol  thia  report) 

Unclassified 

IS*.  OECL  ASSIFICATION  'DOWNGRADING 
SCHEDULE 

*6  DISTRIBUTION  STATEMENT  (ol  Ihl a  Rtport) 


Approved  for  public  release;  distribution  unlimited. 


17  DISTRIBUTION  STATEMENT  (ol  the  abetroct  entered  In  Block  20.  II  dlllerent  Iron*  Roport) 


is  supplementary  notes 


IS  KEY  WOROS  (Continue  on  reverse  aide  it  nec  eaamey  and  Identity  by  block  number) 

Hygrothermal  effects  Temperature  Composite  materials 

Moisture  content  Stress  Analysis  Transient  response 

Elastic  properties  Finite  element  analysis  Hydroelasticity 


20  ABSTRACT  (Continue  on  reverse  atdm  If  net  misery  ond  Identify  by  block  number) 


(SEE  REVERSE  SIDE) 


1— - - -  I 

DO  *  JAN  71  1473  EOlTION  O*  I  NOV  AS  (S  OBSOL*  Tfc 

UNLLASS I  FI LD 


security  Classification  of  rMt$  page  <When  Par,  Fmereo 


_ UNCLASSIFIED _ 

MCUWITV  CLAMirtCATIQH  OF  THIS  PAOtfWTin  Datm  Knfrl) _ _ 

Block  No.  20 

ABSTRACT 

A  system  of  coupled  differential  equations  involving  moisture,  tempera¬ 
ture  and  stress  are  derived  within  the  framework  of  classical  thermodynamics. 
Analyzed  are  the  transient  character  of  the  stress  and  energy  density  distri¬ 
bution  around  a  circular  cavity  in  a  cross-ply  graphite  T300/5208  epoxy  resin 
laminate  assumed  to  behave  in  a  quasi-isotropic  manner.  The  laminate  is 
stretched  uniaxially  such  that  the  mutual  interaction  of  moisture,  temperature 
and  stress  depend  on  whether  the  laminate  is  in  a  state  of  plane  strain  or 
plane  stress.  Unlike  the  case  of  moisture-temperature  and  stress  being  uncou¬ 
pled,  the  diffusion  coefficients  for  plane  strain  must  be  distinguished  from 
those  for  plane  stress  when  stress  coupling  effects  are  included.  Time- 
dependent  numerical  results  are  obtained  from  a  finite  element  formulation 
and  reveal  that  stress  coupling  can  have  a  marked  influence  on  the  solution, 
both  qualitatively  and  quantitatively.  This  is  particularly  significant  near 
defects  or  cavities  where  the  stress  and  energy  state  experience  high  elevation. 
The  neglect  of  stress  coupling  can  lead  to  unconservative  predictions.  Pos¬ 
sible  failure  sites  are  predicted  from  the  strain  energy  density  criterion 
which  assume  that  excessive  dilatation  is  associated  with  material  separation 
while  excessive  distortion  with  yielding. Aas  the  magnitude  of  the  applied 
uniaxial  stress  is  increased,  the  predicted  failure  site  tends  to  move  closer 
to  the  cavity  boundary.  The  present  theory  and  method  of  analysis  provide  a 
quantitative  assessment  of  moisture,  temperature  and  mechanical  load  effects 
<?n  material  failure. 


_ UNCLASSIFIED _ 

JCCUrtfTv  CLASSIFICATION  of  This  PAGE'Whtn  Dtit  Emtrtd: 


Foreword 


The  work  In  this  report  is  prepared  for  the  Army  Materials  and  Mechanics 
Research  Center  (AMMRC)  under  Contract  No.  DAAG46-82-K-0006  with  the  Institute 
of  Fracture  and  Solid  Mechanics  at  Lehigh  University,  Bethlehem,  Pennsylvania. 
The  authors  are  particularly  grateful  to  the  Project  Manager,  Mr.  J.  F.  Dignam 
and  Technical  Monitor,  Or.  S.-C.  Chou  from  AMMRC.  Their  constant  interest  and 
encouragement  have  made  the  progress  in  this  work  possible. 


Accession  For 

KITS  GRAM 
DTI C  TAB 
Unannounced  F] 

Justification.— 


|  By - - 

^Distribution/ 

Availability  Coio3 
iAvail  and/or 
Dlat  |  Special 


A  GENERALIZED  COUPLE  THEORY  OF  HYGROTHERMAL-ELASTICITY: 
TRANSIENT  EFFECTS  IN  COMPOSITE  LAMINATE  WITH  CIRCULAR  CAVITY 


by 

G.  C.  Sih  and  M.  T.  Shih 
Institute  of  Fracture  and  Solid  Mechanics 
Lehigh  University 
Bethlehem,  Pennsylvania  18015 

ABSTRACT 

A  system  of  coupled  differential  equations  involving  moisture,  temperature 
and  stress  are  derived  within  the  framework  of  classical  thermodynamics.  Ana¬ 
lyzed  are  the  transient  character  of  the  stress  and  energy  density  distribution 
around  a  circular  cavity  in  a  cross-ply  graphite  T300/5208  epoxy  resin  laminate 
assumed  to  behave  in  a  quasi -isotropic  manner.  The  laminate  is  stretched  uni¬ 
axial  ly  such  that  the  mutual  interaction  of  moisture,  temperature  and  stress  de¬ 
pend  on  whether  the  laminate  is  in  a  state  of  plane  strain  or  plane  stress.  Un¬ 
like  the  case  of  moisture- temperature  and  stress  being  uncoupled,  the  diffusion 
coefficients  for  plane  strain  must  be  distinguished  from  those  for  plane  stress 
when  stress  coupling  effects  are  included.  Time-dependent  numerical  results  are 
obtained  from  a  finite  element  formulation  and  reveal  that  stress  coupling  can 
have  a  marked  influence  on  the  solution,  both  qualitatively  and  quantitatively. 
This  is  particularly  significant  near  defects  or  cavities  where  the  stress  and 
energy  state  experience  high  elevation.  The  neglect  of  stress  coupling  can  lead 
to  unconservative  predictions.  Possible  failure  sites  are  predicted  from  the 
strain  energy  density  criterion  which  assume  that  excessive  dilatation  is  asso¬ 
ciated  with  material  separation  while  excessive  distortion  with  yielding.  As  the 
magnitude  of  the  applied  uniaxial  stress  is  increased,  the  predicted  failure  site 
tends  to  move  closer  to  the  cavity  boundary.  The  present  theory  and  method  of 


analysis  provide  a  quantitative  assessment  of  moisture,  temperature  and  mechani¬ 
cal  load  effects  on  material  failure. 

INTRODUCTION 

A  critical  aspect  of  composite  materials  is  the  influence  of  environments 
such  as  temperature  and  moisture  in  addition  to  mechanical  loads  on  their  life 
expectancy.  Despite  the  numerous  test  data  collected  in  the  past,  the  degrada¬ 
tion  of  composite  stiffness  and  strength  due  to  environment  is  still  not  under¬ 
stood  and  has  not  been  appropriately  quantified.  Empirical  approaches  are  not 
only  costly  but  cannot  lead  to  a  fundamental  understanding  of  coupling  between 
mechanical  deformation  and  diffusion.  In  the  absence  of  this  knowledge,  compos¬ 
ites  cannot  be  successfully  applied  in  design  and  manufacture. 

In  a  series  of  papers  [1-4],  coupled  diffusion  theories  involving  moisture 
and  temperature  were  derived  on  the  basis  of  classical  thermodynamics.  Numerous 
examples  were  solved  to  show  that  uncoupling  of  moisture  and  temperature  cannot 
be  justified  in  general.  Appreciable  difference  attributed  to  coupling  can  re¬ 
sult  when  the  surface  temperature  is  changed  suddenly.  This  alters  the  transient 
character  of  the  moisture  distribution  that  also  affects  the  resulting  stresses 
and/or  strains.  Loading  type  also  interacts  with  the  diffusion  process  in  a  com¬ 
plex  fashion.  Because  of  the  complexities  of  the  theoretical  treatment  and  pro¬ 
cedures  involved  in  determining  the  physical  constants,  only  a  few  cases  have 
been  analyzed  up  to  this  date. 

A  recent  experiment  [5]  has  shown  that  the  applied  stress  can  indeed  signifi¬ 
cantly  affect  the  moisture  and  temperature  distribution  in  composite  materials. 
Such  findings  are  not  unexpected  and  suggest  the  need  to  consider  coupling  between 
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mechanical  deformation  and  diffusion  by  moisture  and  temperature.  Such  a  theory 
leading  to  three  simultaneous  partial  differential  equations  will  be  presented 
in  this  work.  The  method  of  analysis  will  involve  an  iteration  procedure  in 
which  initial  values  of  moisture  and  temperature  will  be  assumed  to  calculate 
the  corresponding  stress  state.  The  iteration  is  repeated  many  times  until  all 
governing  equations  and  transient  boundary  conditions  are  satisfied.  Results 
are  obtained  specifically  for  a  T300/5208  graphite  epoxy  laminate  with  a  circu¬ 
lar  hole  stretched  uniaxially  while  the  surface  moisture  is  changed  suddenly  at 

a  constant  surface  temperature.  The  stresses  near  the  hole  are  found  to  oscil- 

★ 

late  in  time  and  to  deviate  significantly  from  the  corresponding  uncoupled  re¬ 
sults  for  small  time.  In  accordance  with  the  strain  energy  density  criterion 
[6,7],  the  uncoupled  solution  can  underestimate  the  threshold  of  material  damage 
due  to  fracture  initiation. 

A  COUPLED  THEORY  OF  HYGROTHERMAL  ELASTICITY 

Based  on  the  classical  thermodynamic  argument,  a  coupled  theory  of  mechani¬ 
cal  deformation  and  diffusion  involving  moisture  and  temperature  will  be  derived. 
The  governing  equations  reduce  to  those  treated  earlier  when  deformation  is  un¬ 
coupled  [1-4].  In  the  absence  of  moisture,  the  classical  coupled  equations  of 
thermoelasticity  are  obtained. 

Thzmo dynamic.  Stctfz.  Let  the  state  of  a  unit  volume  of  a  solid  be  described 
by  the  temperature  T,  moisture  C,  entropy  S,  chemical  potential  P,  stress  tensor 
a  and  strain  tensor  e.  The  initial  state  of  the  solid  is  assumed  to  be  in  thermo¬ 
dynamic  equilibrium  with  the  properties  TQ,  CQ,  SQ,  PQ  and  vanishing  g  and  e. 

★ 

In  this  work,  uncoupling  refers  only  to  the  interdependence  of  mechanical  defor¬ 
mation  and  diffusion.  Moisture  and  temperature  are  assumed  to  be  coupled. 
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For  small  deviations  from  the  equilibrium  state,  a  potential  $  may  be  expanded 
in  a  power  series  by  neglecting  those  terms  higher  than  the  second  order  [8]: 


VT  +  Vc  -3'ii 


1 
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in  which  p  is  the  mass  density,  cp  the  specific  heat  capacity,  a  the  linear  ther¬ 
mal  expansion  coefficient,  8  the  linear  moisture  expansion  coefficient,  G  the 
shear  modulus  of  elasticity  and  k  =  E/3(l-2vp)  the  bulk  modulus  of  elasticity 
with  E  being  the  Young's  modulus  and  vp  the  Poisson's  ratio.  In  equation  (1), 
the  coefficient  dc  and  dt  are  defined  as 


-  f  “^1 
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The  chemical  potential  p  is  given  by 


(2) 


p  "  f  "  Po  +  dtAT  +  dcAC  -  6Gii 


(3) 


The  entropy  S  and  strain  components  e •  •  are  obtainable  from  the  thermodynamic 

*  J 

potential  $  according  to 


S  =  ‘  If  =  So  +  AT  •  dtAC  +  a0ii 


(4) 


and 
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The  coupling  of  moisture,  temperature  and  mechanical  deformation  is  thus  reflected 
in  $  as  given  in  equation  (1). 


CorueAvation  Lam.  The  conservation  of  mass  and  energy  will  be  invoked  to 
describe  the  diffusion  and  stress  state  in  the  solid  at  a  given  time  t.  Let  f 
be  the  flux  of  moisture  per  unit  area  and  time,  mass  conservation  requires  that 

vf|§=0  (6) 


Moreover,  the  conservation  of  energy  leads  to  an  expression  for  the  rate  of  change 
of  internal  energy  with  time: 
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in  which  q  stands  for  the  heat  flux  vector  and 


v  .  q  = 


TH* iP 


3C 

at 


(8) 


where  AP  =  P-PQ.  Application  of  the  Onsager  principle  [2,3]  leads  to  the  follow¬ 
ing  expressions  for  the  moisture  and  heat  flux  vectors 


with  L-.  (i,j  =  1,2)  being  the  Onsager  coefficients  such  that 

'  J 

I~1 1  >  0,  1*22  >  0*  L'j 2  =  l-2i  <  1  ^22 

The  vectors  and  are  given  by 

xf  =  -  TV  (y)  •-  (dtAT  +  dcAC  -  0o11)  y^  -  dtvT  -  dcvc  +  6Vaii 

?,  *  •  T  7T 

Substituting  equations  (11)  into  (9)  and  making  use  of  equation  (6)  yield 

(Lndt  +  h2»v2T  +  LnV2c  -  Lu672oii  ■  £ 

Similarly,  combining  equations  (4),  (8)  and  (9)  renders  the  relation 


<L21dt  *  T-  L22>72T  +  L21dc72C  -  L21672°1i 
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In  deriving  equations  (12)  and  (13),  it  was  assumed  that  |dcAC|  and  lea^l  ar 
smaller  in  comparison  with  | d tTQ |  and  that  |aT|  «  TQ.  The  remaining  coupled 
equation  will  be  found  from  the  theory  of  elasticity. 


Equation  Motion.  The  equation  of  motion  can  be  written  as 


where  u.  are  the  components  of  the  displacement  vector  u.  From  the  linear  the¬ 
ory  of  elasticity,  th„  stress-strain  relation  is 


v  1  +v 

0ij  =  2G  [eij  +  T^T  ekk5ij  "  7^  (aAT  +  6AC)«iJ] 


where  vp  is  the  Poisson's  ratio  of  the  solid.  Putting  equation  (15)  into  (14) 
and  knowing 


en  =  7  •  s 


it  can  be  shown  that 


1 "V  32a,, 
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2  p  2  p  p 


For  a  quasi-static  process  where  the  speed  of  diffusion  and  load  application  is 
much  slower  than  the  speed  of  shear  wave  velocity 


c? .  /QL)  a 
2  '1-Zv  '  P 


equation  (17)  may  be  simplified  as 


1  -v 

V2ain.  +  aV2T  +  8V2C  =  0 


Solving  for  from  equation  (19),  it  is  obvious  that 


0  =  .  [aAT  +  8AC]  +  4 ) 


in  which  the  Laplacian  of  <j»  vanishes,  i.e.. 


v24>  =  0 


It  is  seen  that  equations  (12),  (13)  and  (19)  are  the  system  of  governing  equa¬ 
tions  solving  for  T,  C  and  . 

Coupled  Equations.  The  system  of  equations  derived  earlier  can  be  rearranged 
into  forms  that  are  more  recognizable  when  compared  to  the  corresponding  equations 
in  the  uncoupled  theories.  Eliminating  in  equations  (12)  and  (13)  by  applying 
equation  (20),  the  following  coupled  equations  are  found: 


V2<J>  =  0 


in  which 


R  =  T  +  N<|> 


Note  that  there  are  five  physical  constants  D,  V,  X,  v  and  N  in  this  theory.  They 
are  interconnected  through  the  parameters  L .  .  (i,j  =  1,2),  d  ,  dt,  etc.,  as  fol- 

I  J  L.  L 

lows: 


D  =  (l-Xv)D 


V  =  (l-Xv)Dk 


in  which  Dm  and  are  given  by 
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The  quantity  A  stands  for 


A  =  +  d 
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The  remaining  constants  in  equations  (22)  and  (23)  are 


<AL„  +  L,2)(fE  -  f^)/(AT0Ln  +  ^  L 22  +  2L12ATo> 
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This  completes  the  formulation  of  the  coupled  theory  of  hygrothermal  elasticity. 


PLANE  THEORIES  OF  HYGROTHERMAL  ELASTICITY  AND  SPECIAL  CASES 


When  a  length  dimension  of  the  solid,  say  in  the  z-direction,  is  either  very 
small  or  large  in  comparison  with  the  planar  dimensions  in  the  xy  domain,  certain 
approximations  on  the  governing  equations  can  be  made.  This  will  not  only  alter 
the  form  of  the  equations  but  also  the  definitions  of  the  physical  constants. 

The  two  cases  that  will  be  discussed  are  plane  stress  and  plane  strain  state  as 
defined  in  the  classical  theory  of  elasticity. 


Plane.  S£n<u>4.  For  thin  plates  with  surfaces  free  from  tractions,  the  trans¬ 
verse  normal  stress  component  033  may  be  assumed  to  vanish  and  o^.  equals  to 
all+cr22  ™  e<luat*Qns  (12)  and  (13).  Since  it  is  easier  to  perform  experiments 
on  thin  plates,  the  physical  constants  D,  P,  D  ,  Dh,  etc.,  in  equations  (24) 
through  (27)  shall  refer  to  the  case  of  plane  stress. 


Plane  Strain.  For  long  cylindrical  bodies  with  axes  along  the  z-di recti  on. 


the  longitudinal  strain  components  may  be  assumed  to  vanish,  i.e.. 
Making  use  of  equation  (15)  with  633  =  0,  it  is  found  that 


a33  =  vp  ^  ai  l  +a22  ^  ~  +  SAC] 

Substituting  equation  (28)  into  (15)  gives 


Here,  akk  =  °yi+°22‘  Equations  (29)  and  (15)  are  equivalent  if 


(28) 


(29) 


a*  =  (l+vp)a,  8*  =  (1+Vp)g 


(30) 


Henceforth,  those  quantities  referred  to  the  plane  strain  state  will  be  distin¬ 
guished  by  attaching  an  asterisk  to  them: 


Thus,  the  same  governing  equations  apply  to  plane  stress  and  plane  strain.  Only 
the  constants  need  to  be  replaced  as  indicated  in  equations  (30)  and  (31). 

Uncoupled  Vz^onmation.  When  the  mechanical  deformation  term  N$  in  equation 
(23)  is  neglected,  equations  (22)  reduce  to 


X  JI 
A  3t 


which  have  already  been  discussed  previously  [2], 


Thzxjnoztcaticlty  Thzony.  By  dropping  those  terms  in  equations  (12),  (13)  and 
(15)  pertaining  to  moisture  change,  the  coupled  thermoelastic  equations  may  be 
obtained: 


“p  S  *  °To  ir 3  <“*1  +  r  42>72t 
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aij  =  26  [eij  +  rar  ekk6ij  '  aAT6ij^ 


The  coefficients  in  equations  (33)  may  be  redefined  to  coincide  with  those  given 
in  the  text  book. 


DETERMINATION  OF  PHYSICAL  CONSTANTS 


One  of  the  major  tasks  of  evaluating  the  coupling  effects  of  moisture,  tem¬ 
perature  and  mechanical  deformation  is  a  knowledge  of  the  physical  constants. 

The  material  depicted  will  be  that  of  a  T300/5208  graphite  epoxy  composite  lami¬ 
nate  which  will  be  assumed  to  behave  as  a  quasi-isotropic  solid.  The  coupled 
theory  solution  will  be  obtained  for  the  case  where  the  coupling  effects  are 
weak  such  that  the  results  can  be  compared  with  the  uncoupled  experimental  data 
to  extract  the  coupling  constants. 

Rule,  oh  UixtuAe..  The  rule  of  mixture  will  be  first  applied  to  obtain  an  aver¬ 
age  value  of  the  thermal  expansion  coefficient  a  for  the  composite  with  56%  fi¬ 
ber  volume  fraction  V*.  Let  a  and  a*  denote  the  respective  values  of  the  thermal 

t  m  r 

expansion  coefficient  for  the  matrix  and  fiber.  The  value  a  for  the  composite 
may  be  estimated  from  [10] 


The  composite  modulus  in  the  direction  of  fiber  reinforcement  may  be  approxi¬ 


mated  by 


’ll 


*  Efvf  + 


EV 
m  m 


(37) 


The  modulus  in  the  transverse  direction  will  be  denoted  by  E22.  An  approximate 
value  E  for  a  laminate  with  cross-plies  can  be  written  as  [12] 


E  = 


_  3 


*11 


5  E 
8  fc22 


=  67.155  GPa 


(38) 


Based  on  these  assumptions,  the  quasi-isotropic  properties  of  the  T300/5208  graph¬ 
ite  epoxy  composite  laminate  are  given  in  Table  1. 


Table  1  -  Material  properties  for  T300/5208  graphite  epoxy  composite  laminate 

_ Cli 1 _ 


Graphite 

Fiber 

Epoxy 

Matrix 

Composite 

Lami nate 

Axial  Modulus 

En(GPa) 

276.00 

E*: 

3.51-0.003T 

0.142C 

156.08 

Transverse  Modulus 

E22(GPa) 

13.8 

E": 

3.51 -0.003T 
0.142C 

13.80 

Poisson's  Axial  Ratio 

<Vn 

Poisson's  Transverse  Ratio 

(V22 

0.20 

0.25 

0.33 

0.33 

— 

|  0.33 

Axial  Thermal  Coefficient 

“ll/°C 

-  1.8  x  10'° 

45.0 

x  10"6 

Transverse  Thermal  Coefficient 
a22/<>C 

21.6  x  10"6 

45.0 

x  10'6  J 

>31.3  x  10"6 

Moisture  Coefficient  - 

B/wt%H20  0  2.68  x  10"J  2.68  x  10" 


The  specific  heat  cp  of  the  T300/5208  laminate  can  be  estimated  from  those 
for  the  matrix,  water  and  fiber.  They  are  given  by 

(cn)  =  0.22  cal/g°C,  (cn)  =  0.99  cal/g°C 
p  m  P  H20 

(cp)  =  0.336  +  9.50  x  10‘5  T  -  1.70  x  104  (39) 

=  0.167  cal/g°C  at  T  =  294°K 

It  is  also  known  from  [14]  that  the  weight  of  one  cubic  meter  of  T300/5208  composite 
is  1,590  Kg  and  1,608.48  Kg  after  fully  saturated  by  moisture.  This  yields  a  dif¬ 
ference  of  18.48  Kg.  Since  the  fiber  weight  percentage  is  73.5,  the  fiber  weight 
percentage  corresponding  to  the  fully  saturated  state  may  be  computed  as 

(%  wt)f  =  T | 606?46  (73'5%)  =  72*7%  (40) 

Under  these  considerations,  a  value  of  cp  for  the  laminate  is  obtained: 

c  =  {%  wt)f(c  )  +  (%  wt)  (c  )  +  (%  wt)H  n(c  ) 
p  f  P  f  m  P  m  H2°  p  H20 

=  (0. 727 ) (0. 1 67 )  +  (0. 262) (0.220)  +  (0.012) (0.99) 

=  0.19041  cal/g°C  (41) 

Unit  ConvoAAion.  Before  deriving  the  physical  constants  in  equations  (22)  and 
(23),  it  is  necessary  to  make  use  of  the  results  in  [3,4]  which  refer  to  C  as  mass 
of  moisture  per  unit  volume  of  void.  The  corresponding  constants  D  ,  V  and  X  v 


4- 


for  the  T300/5208  material  are  given  by 


a  0.1 ,  Xyvy  =  0.25 


The  quantities  D,  P,  and  Av  in  equations  (22)  correspond  to  C  expressed  as  mass 
of  moisture  per  unit  volume  of  solid.  The  relationships 


1  +  ^  Vv 


'"Vv 


_  W  V 

1  +  IT  Vv 

given  in  [2]  may  be  used  for  the  conversion.  From  equations  (42)  and  (43),  it  is 
found  that 


$  *  0.07317,  Av  =  0.02439 


The  individual  values  of  D,  P,  A  and  v  may  be  obtained  by  applying  equations  (24) 
and  the  relations 


D  C  Q  Q 

,  _  m  oyo  yo 

K  W  v ' 

n  g  o  p 


in  which  p  is  the  mass  density,  Rg  the  gas  constant,  CQ  the  weight  of  moisture  ab¬ 
sorbed  by  one  cubic  meter  of  the  composite  at  21 °C  and  75%  RH,  Qq  is  a  constant, 
Tq  *  294°K,  and  Cp  is  given  by  equation  (41). 


Coupled  Constant*  fax  Plane.  Stress.  It  was  pointed  out  in  [2]  that  the  dif¬ 
fusion  coefficients  for  sudden  moisture  change  and  should  be  distinguished 

from  and  corresponding  to  sudden  temperature  change.  These  quantities  are 

★ 

expressible  in  terms  of  D,  V,  A  and  v  as  follows  : 

Dc  a  l-Xvf(u.Av)  D 

m  1-Av 

(46) 

Dc  =  0 

h  l-Av 


and 


0T  .  f(utpl  v 

m  l-Av 

(47) 

0T  ,  l->vf|Vu,,XvI  „ 

in  which  the  function  f(u,Av)  stands  for 

f(u,Av)  =  [1  +  u  +  v'4u(l-Av)  J-1  =  0.62207  (48) 


since  the  numerical  values  of  u  =  D/P  and  Av  are  known  from  equations  (44).  Know¬ 
ing  from  [14]  that 


o£  *  D„  exp ( E  / R  T  )  =  7.7954  x  10"7  cm2/hr 
m  o  o  g  o 


(49) 


As  discussed  in  [2],  these  constants  can  be  measured  experimentally  by  subjecting 
the  specimen  to  sudden  moisture  or  temperature  change. 
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where  EQ  is  the  activation  energy.  The  moisture  diffusivity  in  equation  (49) 
for  the  T300/5208  material  is  calculated  from  the  following  constants: 


DQ  *  1.53  x  103  cm2/hr,  EQ  =  1.25  x  104  cal/g  •  mole 
R  =  1.987  cal/g  .  mole  •  °k,  T  =  294°K 


(50 


Inserting  the  value  of  Dm  into  the  first  of  equations  (46)  yields  D  and  together 
with  the  first  of  equations  (44)  gives  P,  i.e., 

D  =  7.72242  x  10'7  cm2/hr,  V  =  1.05541  x  10"5  cm2/hr  (51) 


There  remains  the  determination  of  A  and  v.  Combining  the  second  of  equations 
(44)  with  (45),  the  constant  Qq  can  be  solved: 


/ 


pcpXvRTo 
uT - 


=  9.55650  x  105  m2/s2 


(52) 


in  which  CQ  =  18.484  Kg  and  p  =  1,590  Kg  for  the  T300/5208  composite  [14].  The 
second  of  equations  (44)  and  (45)  now  yield 


„-3 


A  =  0.03235  kg  .  m”7°C,  v  =  0.75389°C/kg 


m 


-3 


(53) 


This  completes  the  numerical  evaluation  of  D,  V,  A  and  v  for  the  case  of  plane 
stress. 


It  follows  immediately  from  equations  (24)  that 


(54) 


Dm  =  7.91547  x  10  '  cm2/hr,  Dh  =  1.08179  x  10"3  cm2/hr 

When  deformation  is  coupled  with  diffusion,  the  additional  constant  N  given  by 
the  third  of  equations  (27)  must  be  evaluated.  Knowing  that  a  =  31.3  x  10"®  and 
Vp  =  0.33  as  given  in  Table  1,  and  E  =  67.155  GPa  and  Cp  as  computed  by  equations 
(38)  and  (41 ),  the  value 

N  =  7.60997  x  10"9  ms2oC/Kg  (55) 

is  evaluated  for  plane  stress.  A  glance  at  the  six  expressions  in  equations  (25) 
to  (27)  show  that  the  remaining  six  unknowns  dc,  dt.  A,  L^,  L12  and  l_22  can  be 
obtained  without  difficulty.  Their  values  are 

d  =  1,230.33  m5/Kgs2 

dt  =  -3,330.931  m2/s2oC  (56) 

A  =  -2,257.93  m2/s2°C 

and  the  Onsager  coefficients  are  given  by 

L-ji  =  3.13576  x  10"18  Kg-s/m3 

L12  =  4.93980  x  10"12  Kg/m-s  (57) 

L22  =  7.19903  x  10"4  Kg-m/s3 


The  above  numerical  evaluation  is  based  on  the  best  approximation  by  assuming  that 


L22  remains  unchanged  as  Dm,  D^,  X  and  v  are  matched  with  those  obtained  from 
the  classical  theory  where  moisture  and  temperature  are  uncoupled.  Alternatively, 


L.J2  may  be  assumed  to  remain  unchanged  which  will  lead  to  a  slightly  different 
value  of  A. 

Coupled  Constant*  hoti  Plane.  S&iain.  The  conversion  to  plane  strain  can  be 
made  in  a  straightforward  manner.  First,  the  constants  a*,  b*,  v*  and  E*  can  be 
computed  from  equations  (30): 


a*  =  4.15353  x  10_5/°C,  B*  =  3.5912  x  10'3/wt%  HgO 
v*  =  0.493,  E*  =  75.9328  x  109  N/m2 

With  the  aid  of  equations  (31),  the  remaining  quantities  are  also  obtained: 


(58) 


cj  =  785.1243  m2/s2°C 

dj  =  7,013.8  m5/Kgs2  (59) 

d*  =  -3,330.6  m2/s2oC 


Following  the  definition  given  by  the  third  of  equations  (27),  N*  is  obtained 


N*  =  a* 


2E*a*2n_1 

T-“J 

p 


=  1.12402  x  10“8  ms2°C/Kg 


(60) 


The  coupling  constants  are  determined  by  matching  the  theoretical  boundary  value 
solution  for  a  problem  with  very  weak  coupling  to  the  experimental  data  without 
the  influence  of  coupling.  Since  the  physical  constants,  once  determined,  are 
independent  of  boundary  conditions,  they  can  be  used  in  general  for  problems  where 
coupling  influence  is  significant. 


Similarly,  equations  (25)  in  terms  of  plane  strain  render 


D*  =  +  d*)Ln  =  2.59555  x  10-6  cm2/hr 

°h  -  <A*Tohl  +  T-  L22  +  2Li2)(A*/T0)/<^£  - 

0  op 

=  1.26108  x  10' 5  cm2/hr 

in  which  A*  is  obtainable  from  equation  (26)  as 

A*  =  2E^”*f*  +  d*  =  -  392.183  m2/s2-°C 
~VP 

Therefore,  X*  and  v*  follow  from  equations  (27): 

X*  =  ( L-j i A*  +  j-  L12)/D*  =  0.04454  Kg  •  m“3/°C 

v*  =  (f-r^  +  A*)/(^£  -  2-E*”*2)  =  1.34389  °C/Kg  •  m'3 
o  L1 1  o  p 

Finally,  D*  and  V*  can  be  calculated  as 
D*  =  (l-X*v*)D*  =  2.44018  x  10"6  cm2/hr 
V*  =  (i _x*v*)D*  =  1.18559  x  10"5  cm2/hr 
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(61 


(6^ 


The  constants  in  equations  (63)  and  (64)  will  be  used  subsequently  for  analyzing 
the  transient  hygrothermal  stresses  around  a  circular  hole  in  a  stretched  plate. 


TRANSIENT  DIFFUSION  AND  STRESS  DISTRIBUTION  NEAR  CIRCULAR  HOLE 


Consider  the  problem  of  a  circular  hole  of  radius,  a,  in  a  square  plate  with 
dimensions  LxL  as  shown  in  Figure  1.  The  boundary  conditions  on  the  hole  con¬ 
sist  of  sudden  rise  in  moisture  or  temperature  and  are  independent  of  the  polar 
angle  0.  A  uniform  stress  of  magnitude  aQ  is  applied  in  the  y-direction  on 
such  that  only  the  shaded  area  need  to  be  analyzed.  It  is  of  interest  to  investi¬ 
gate  the  interaction  between  the  stresses  induced  by  diffusion  and  those  by  me¬ 
chanical  loading. 

Finite.  Element  Analy^i*.  The  basic  time-dependent  finite  element  formulation 
follows  that  in  [15]  for  the  case  when  the  deformation  is  uncoupled  from  the  dif¬ 
fusion  process.  The  governing  differential  equations  correspond  to  those  in 
equations  (32)  which  are  similar  to  equations  (22)  except  for  R  being  replaced  by 
T  as  <f>  =  0.  Without  going  into  details,  it  should  be  mentioned  that  the  varia¬ 
tional  calculus  formulation  in  [15]  is  unique  because  the  time  dependent  portion 
of  the  problem  is  solved  analytically  as  an  eigenvalue  problem  by  application  of 
Laplace  transform  in  the  time  variable.  This  choice  is  also  adopted  in  this  work 
and  will  not  be  elaborated  furtner.  Referring  to  Figure  2  for  the  finite  element 
grid  pattern  layout  in  the  first  quadrant.  A  total  of  139  triangular  elements 
and  87  nodes  are  used.  Computations  are  carried  out  by  letting  a  =  1.0  and  L 
=  4.0. 


For  the  convenience  of  numerical  computation,  the  first  two  expressions  of 
equations  (22)  are  cast  into  dimensionless  forms: 


v2C  =  (C-x*R) 

e  3t 


—  V2R  =  |-  (R-v*C) 

£  3t 


(65) 
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Figure  1  -  Circular  hole  in  a  stretched  plate 


in  which  x  =  et/h2  is  the  normalized  time  parameter  with  e  being  a  time  scale 
factor  and  h  the  plate  thickness.  The  space  variables  in  the  Laplacian  operator 
v2  are  normalized  with  reference  to  h,  i.e.,  x/h,  y/h  and  z/h.  The  time  scale 
factor  e  is  chosen,  respectively,  as  10“^  cm2/hr  and  1.0  cm2/hr  for  suddenly  ap¬ 
plied  moisture  and  temperature  boundary  conditions. 


Suddzn  MoiAtuAz  ZUz.  Let  the  initial  conditions  of  moisture  and  temperature 
on  Tj  and  be  such  that 

ACg ,  on  r j 


0  ,  on  r. 


AT  =  0,  on  rT  and  rTT 


More  specifically,  the  moisture  and  temperature  on  are  kept  constant  at  0% 

RH  and  21°C,  respectively,  at  all  time  t  while  the  moisture  on  Tj  is  raised  from 
0%  to  75%  RH  at  constant  temperature.  Results  for  a  constant  applied  stress  of 
oQ  *  5  x  10^  N/m2  are  first  presented  while  the  influence  of  raising  oQ  to  10® 
N/m2  will  be  discussed  only  in  terms  of  the  strain  energy  density  function. 

Let  the  average  moisture  and  temperature  be  weighed  through  the  volume  of 
each  triangular  element  as  follows: 


C  =  1  /  Cdv,  T  -  l  /  Tdv  (67) 

In  the  sequel,  all  solid  curves  refer  to  stress  coupling  while  the  dotted  curves 
correspond  to  stress  uncoupling.  Moisture  and  temperature  are  always  coupled. 
Figure  3  gives  a  plot  of  (£-C0)%  dry  weight  against  t.  The  coupled  solution  can 


0.0  2.0  4.0  6.0  8.0  10. 


Figure  3  -  Average  moisture  gain  for  sudden  moisture  change  on 
circular  cavity  from  21 °C  0%  to  75%  RH.  Material 
T300/5208  graphite/epoxy  laminates 


be  40%  to  50%  higher  than  the  uncoupled  solution.  The  average  temperature  change 
(T-T0)  as  a  function  of  t  is  displayed  in  Figure  4.  Note  that  (T-T  )  decreases 
at  first  and  then  rises  sharply  while  the  corresponding  results  without  stress 
coupling  remains  virtually  constant  with  respect  to  time.  Figures  5  and  6  ex- 


aQ  =  5  x  107  N/m2 


-  coupled 

-  uncoupled 


Figure  4  -  Average  temperature  rise  for  sudden  moisture  change  on 
circular  cavity  from  21 °C  0%  to  75%  RH.  Material 
T300/5208  graphite/epoxy  laminates 


hibit  the  change  of  (TT-C0)  %  dry  weight  and  (T-T  )  with  the  normalized  radial  dis¬ 
tance  r/a  for  t  =  1.0.  It  is  seen  from  Figure  5  that  the  average  moisture  content 
decreases  with  increasing  distance  r/a.  The  stress  coupling  effect  only  changed 
the  results  quantitatively  and  not  qualitatively  in  that  the  solid  curves  are  al¬ 
ways  higher  than  the  dotted  curves  for  a  given  r/a  value  but  the  remained  constant 
for  all  0.  The  uniaxial  stress  level  of  aQ  =  5  x  107  N/m2  is  not  strong  enough 
to  influence  the  asixymmetric  character  of  C-CQ  which  is  dominated  by  the  initial 
moisture  and  temperature  boundary  conditions.  On  the  other  hand,  the  average 

temperature  change  (T-Tq)  in  Figure  6  is  seen  to  be  strongly  affected  by  stress 
★ 

This  transient  effect  becomes  more  pronounced  when  the  temperature  difference 
(T-T  )  at  the  nodal  points  are  plotted  as  a  function  of  time  [16]. 


0.0  30.  0  (degree) 

Figure  5  -  Moisture  distribution  versus  radial  distance  and 
angle  at  t  =  1.0  for  sudden  moisture  change  on 
circular  cavity 

coupling  both  quantitatively  and  qualitatively.  Axisymmetry  of  the  initially 
prescribed  temperature  conditions  is  clearly  disturbed  with  (T-T  )  exhibiting 
minima  between  e  =  50°  and  60°  as  r/a  is  varied  from  1.36  to  1.66.  The  varia¬ 
tions  of  moisture  and  temperature  with  distance* r/a  along  the  x-axis  are  shown 
in  Figures  7  and  8.  While  the  moisture  change  in  Figure  7  decreases  monotonical- 
ly  with  r/a  with  increasing  deviation  between  the  coupled  and  uncoupled  solution 
as  t  increased,  the  temperature  change  in  Figure  8  tends  to  oscillate  reaching  a 
maximum  near  the  hole  and  then  decay  rapidly  to  the  prescribed  value  on  the  bound 


ary  rn. 

★ 

The  distance  r  and  angle  9  are  average  over  two  adjacent  triangular  elements  en- 
closed  by  the  radial  and  circular  lines  in  Figure  2. 


0.0  30.  0  (degree)  60.  90. 

Figure  6  -  Temperature  distribution  versus  radial  distance  and 
angle  at  t  *  1.0  for  sudden  moisture  change  on 
circular  cavity 

The  time-dependent  stress  components  ax>  ay,  az  and  cxz  are  plotted  as  a 
function  of  r/a  for  different  time  t  in  Figures  9  to  12  inclusive.  All  the  nor¬ 
mal  stresses  a  ,  a  and  a  can  be  tensile  or  compressive  depending  on  the  time 
X  y  z 

elapsed  and  they  do  exhibit  a  significant  difference  from  the  uncoupled  theory. 
The  stress  component  cx  in  Figure  9  is  seen  to  be  predominantly  compressive  for 
x  >  0.24.  It  became  tensile  for  t  =  0.04  at  r/a  =  1.75  and  diminished  to  zero  on 
Tjj.  The  maximum  value  of  ax  at  t  =  7.5  is  -42.5  N/m2  for  the  coupled  case  and 
-37.5  N/m2  for  the  uncoupled  case.  Referring  to  Figure  10,  oy  is  mostly  tensile 
except  for  those  locations  near  the  hole  boundary.  The  applied  mechanical  stress 
aQ  *  5  x  107  N/m2  is  seen  to  have  little  influence  on  cy  which  rises  to  a  maximum 
near  the  hole  rather  than  decaying  monotonically  with  (^y)max  at  r/a  =  1.0  if  oQ 
were  the  only  external  loading. 


Figure  8  -  Variation  of  temperature  versus  distance  along  x  axis  for 
sudden  moisture  change  on  circular  cavity.  Material 
T300/5208  laminates 
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Figure  9  -  Variation  of  stress  ox  versus  distance  along  x  axis  for 

sudden  moisture  change  on  circular  cavity.  Material 
T300/5208  laminates 


(MN/m2) 


The  behavior  of  az  is  given  in  Figure  11  which  is  highly  compressive  for  small 
values  of  r/a  and  becomes  tensilb  as  r/a  is  increased.  The  shear  stress  compo¬ 
nent  a  in  Figure  12  remained  constant  for  most  of  the  part  of  the  x-axis.  Only 
xy 

small  fluctuations  are  observed  near  the  hole.  In  general,  deformation  coupling 
tends  to  increase  the  magnitude  of  the  local  stresses  depending,  of  course,  on 
the  time  x. 


r/a 


Figure  12  -  Variation  of  shearing  stress  versus  distance  along  x  axis 
for  sudden  moisture  change  on  circular  cavity.  Material 
T300/5208  laminates 


Sudd&n  TzmpeAcutuAz  Rlsz.  Suppose  that  the  temperature  on  the  hole  Tj  is  now 
increased  suddenly  while  the  moisture  is  kept  constant.  The  following  conditions 
then  prevail: 


AC  =  0,  on  Fj  and  r j j 


AT  =. 


ATg,  on  Tj 
0  ,  on  r 


II 


(68) 


Again,  the  moisture,  temperature  and  stress  conditions  on  Tjj  are  maintained  un¬ 
changed  as  initially  specified  at  21°C  and  75%  RH.  The  applied  stress  oq  in  the 

O 

y-di recti on  is  10  N/m2. 

As  the  temperature  on  Tj  is  raised  from  21 °C  to  61 °C  with  the  moisture  kept 
constant  at  75%  RH.  hygrothermal  stresses  will  be  induced  and  fluctuate  with 
time.  The  same  T300/5208  composite  material  properties  as  those  given  earlier 
will  be  used  to  calculate  the  numerical  results. 

Figure  1 3  gives  a  plot  of  the  average  moisture  'hange  (C-CQ)  expressed  in  terms 
of  %  of  dry  weight  of  the  laminate  against  the  dimensionless  time  parameter  i. 

The  time  scale  factor  e  in  x  is  chosen  as  1.0  cm2/hr  fsr  the  present  problem  of 
sudden  temperature  change  as  specified  by  the  conditions  in  equations  (68).  The 
solid  curve  indicates  that  more  moisture  is  first  gained  by  the  composite  when 
stress  coupling  effect  is  taken  into  account.  The  moisture  content  then  decreases 
as  time  increases.  The  peak  of  the  solid  curve  is  higher  than  that  of  the  dotted 
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curve  and  it  occurs  at  a  smaller  value  of  x.  The  average  temperature  (T-T  )  in 


Figure  13  -  Average  moisture  gain  versus  time  for  sudden  temperature 
change  on  circular  cavity  from  21°C  to  61°C  at  75%  RH. 
Material  T300/5208  graphite/epoxy  laminate 


Figure  14  is  seen  to  rise  quickly  to  a  steady  state.  This  occurs  at  x  =  0.25 
and  0.50  for  the  cases  without  and  with  stress  coupling,  respectively.  Similar 


Figure  14  -  Average  temperature  rise  versus  time  for  sudden  temperature 
I  change  on  circular  cavity  from  21°C  to  61°C  at  75%  RH. 

Material  T300/5208  graphite/epoxy  laminate 


plots  are  displayed  in  Figures  15  and  16  to  show  the  variations  of  ( C-C  )  and 
(T-T  )  with  the  polar  angle  0.  Deviations  from  axi symmetry  become  more  pro- 


nounced  for  ( C-CQ )  as  r/a  is  increased  while  the  opposite  trend  is  observed  for 
(T-T0).  It  is  of  interest  to  note  that  (C-CQ)  increases  rapidly  to  a  peak  for 
small  t  near  the  hole,  i.e.,  r  «  a  and  then  it  decreases  sharply  to  zero.  Fig¬ 
ure  17.  This  effect  dies  out  with  time  as  it  is  to  be  expected.  The  average 
temperature  change  (T-T  )  varied  monotonically  with  r/a  for  small  time  and  then 
it  acquires  a  small  peak  for  x  >  0.1  before  it  decays  rapidly.  This  is  shown 
in  Figure  18. 


1.0  2.0  3.0  4.0 

r/a 


Figure  17  -  Variation  of  moisture  versus  distance  along  x  axis 
for  sudden  temperature  change  on  circular  cavity. 
Material  T300/5208  laminate 


The  results  in  Figures  19  to  22  inclusive  pertain  to  plots  of  o  ,  a  ,  a 

a  y  z 

and  a  against  r/a  for  different  time  t.  The  stress  component  a  acquires  a 
xy  x 

peak  near  the  hole  before  it  decreases  in  magnitude  as  indicated  in  Figure  19. 
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Figure  18  -  Variation  of  temperature  versus  distance  along  x  axis  for 
sudden  temperature  change  on  circular  cavity.  Material 
T300/5208  laminate 


Figure  19  -  Variation  of  stress  ax  versus  distance  along  x 

axis  for  sudden  temperature  change  on  circular 
cavity.  Material  T300/5208  laminate 


p 

For  an  applied  stress  level  of  crQ  =  10  N/m2,  stress  coupling  tends  to  lower 

the  peak  value  of  a  .  Figure  20  shows  that  a  decreases  with  r/a  and  the  solid 

x  y 

curves  become  larger  in  magnitude  than  the  dotted  curves  at  r/a  =  1.75.  The 
dominance  of  mechanical  loading  is  clearly  exhibited  when  aQ  is  increased  from 
5  x  10^  N/m2  in  Figure  10  to  10®  N/m2  in  Figure  20.  Note  that  the  maximum  of  a 
now  occurs  at  r/a  =  1.0  and  decays  as  r/a  is  increased.  The  stress  component  a2 
displayed  graphically  in  Figure  21  reveals  that  it  is  compressive  near  the  hole 
and  then  becomes  tensile  at  r/a  >  2.3.  Stress  coupling  does  not  seem  to  influ¬ 
ence  the  shear  stress  axy  as  the  dotted  and  solid  curves  in  Figure  22  almost  co¬ 
incide  with  each  other. 

FAILURE  SITE  PREDICTION:  STRAIN  ENERGY  DENSITY  CRITERION 

The  strain  energy  density  criterion  [6,7]  has  been  used  in  many  applications 
for  predicting  the  sites  of  failure  initiation.  This  is  accomplished  by  comput¬ 
ing  for  the  stationary  values  of  dW/dV  given  by 
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Figure  22  -  Variation  of  stress  a  versus  distance  along  x  axis  for 

xy 

sudden  temperature  change  on  circular  cavity.  Material 
T300/5208  laminates 


dW  _  y 


°ijdeij 


+  f(AT),  i ,j  =  1,2,3 


(69) 


in  which  and  e^-  are  the  rectangular  components  of  the  stress  and  strain 
tensor.  The  function  f(AT)  in  equation  (69)  cannot  be  arbitrarily  dropped  for 
non-uniform  states  of  temperature.  It  is  given  by 


f(4T)  -  Trlr  «2(t-to>2  (70) 

P 

Although  the  function  f (AT)  does  not  affect  the  stress  and  strain  state  in  the 
solid,  it  does  influence  the  energy  density  state  and  hence  the  prediction  of 
failure  sites  should  dW/dV  be  used  in  the  failure  criterion. 

The  basic  idea  in  the  strain  energy  density  criterion  relies  on  the  assump¬ 
tion  that  failure  modes  can  be  uniquely  identified  with  threshold  or  critical 
values  of  dW/dV.  The  locations  of  yielding  and  fracture  initiation  are  asso- 
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dated,  respectively,  with  the  relative  maximum  and  minimum  values  of  dW/dV  with 
reference  to  the  space  variables,  say  r  and  e  in  polar  coordinates  for  two-dimen¬ 
sional  problems.  The  variations  of  dW/dV  with  0  or  r  determine  the  locations  of 
(dW/dV)^^  and  (dW/dV)m^n.  A  critical  value  of  dW/dV  or  (dW/dV)c,  which  can  be 
obtained  from  the  area  under  the  true  stress-strain  curve  in  a  uniaxial  tensile 
test  at  fracture,  is  assumed  to  govern  the  condition  of  fracture  initiation, 
i.e.,  when  (dW/dV)min  =  (dW/dV )c- 

For  a  linear  elastic,  isotropic  and  homogeneous  material,  equation  (69)  can 
be  written  simply  as 


provided  that 

f(AT)  »  j (cxAT)2  (72) 

In  addition,  if  the  problem  is  in  a  state  of  plane  strain,  the  additional  condi 
tion 

°z  =  V°x+V  "  aEAT  (73) 

must  be  imposed.  Hence,  dW/dV  for  the  T300/5208  composite  laminate  can  be  com¬ 
puted  without  difficulty.  Because  of  symmetry  with  respect  to  the  x-axis, 
(dW/dV)m^n  will  occur  at  0  =  0  and  hence  it  is  only  necessary  to  differentiate 
(dW/dV)  with  reference  to  r  for  determining  (dW/dV )mi-n. 


Mo*u£u/lz  Boundary  Condition.  Referring  to  the  boundary  conditions  specified 
by  equations  (66)  where  the  moisture  on  the  hole  is  raised  from  0%  to  75%  RH 
with  the  temperature  held  constant  at  21  °C,  the  variations  of  dW/dV  with  r/a  for 
6*0  are  shown  in  Figure  23  for  t  =  7.5  and  aQ  =  5  x  10^,  10®  N/m2.  Indeed, 


Figure  23  -  Variations  of  strain  energy  density  with  distance  for 
sudden  moisture  change  in  T300/5208  containing  a  circu¬ 
lar  cavity  for  t  =  7.5 


the  solid  curves  possessed  relative  minima  at  r/a  approximately  equal  to  2.975 
7  8 

for  *  5  x  10  N/m2  and  1.850  for  o„  =  10  N/m2.  Increase  in  the  magnitude  of 


the  applied  stress  aQ  tends  to  shift  the  location  of  (dW/dV )ml- n  closer  to  the 
hole  boundary.  Figure  24  gives  a  plot  of  the  minima  of  dW/dV  against  the  dimen¬ 
sionless  time  parameter  t.  For  cQ  =  5.0  x  107  N/m2,  (dW/dV)min  rises  sharply  at 


Figure  24  -  Variations  of  minimum  strain  energy  with  x  for  sudden 
moisture  change  in  T300/5208  containing  a  circular 
cavi ty 

first  and  then  levels  off.  The  time  at  which  failure  occurs  depends  on  the  in¬ 
tersection  of  the  (dW/dV)c  line  with  the  (dW/dV)min  curve.  The  opposite  trend 
is  observed  when  aQ  is  increased  to  10®  N/m2  where  (dW/dV )m^ n  rised  slightly  at 
first  and  then  decayed  monotonically.  When  the  deformation  is  decoupled  from 
diffusion,  the  amplitude  of  (dW/dV)ml.n  decreases  appreciably  as  the  dotted  curve 
in  Figure  25  for  aQ  =  5.0  x  107  M/m2  is  seen  to  lie  below  the  solid  curve. 


Figure  25  -  Minimum  strain  energy  density  as  a  function  of  time  for 
sudden  moisture  change  on  circular  cavity.  Material 
T300/5208  laminates 

TejnpeAAtune.  BoundaA.y  Condition.  The  boundary  conditions  corresponding  to 
sudden  temperature  rise  are  given  in  equations  (68).  A  40°C  of  temperature  in¬ 
crease  from  21 °C  to  61 °C  is  applied  to  the  hole  while  keeping  the  moisture  con¬ 
stant  at  75%  RH.  The  numerical  results  are  displayed  graphically  in  Figures  26 

O 

and  27  for  oQ  *  10  N/m  .  The  general  trend  of  dW/dV  decaying  with  r/a  is  seen 
in  Figure  26  for  x  =  0.2.  A  relatively  weak  minimum  of  dW/dV  is  detected  at 
r/a  =  2.80  and  the  difference  between  stress  coupling  and  uncoupling  is  small. 
Figure  27  gives  a  plot  of  ( dW/dV )ml-n  versus  the  dimensionless  time  parameter  x. 
Aside  from  the  difference  in  magnitude  of  (dW/dV)min  between  the  solid  and  dotted 
curve,  their  trend  is  opposite.  The  coupled  theory  yields  values  of  (dV//dV )mi- n 
increase  with  time  while  the  uncoupled  theory  predicts  the  opposite.  Stress 
coupling  effect  is  thus  seen  to  be  significant  both  quantitatively  and  quali¬ 
tatively. 
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Figure  27  -  Variations  of  minimum  strain  energy  with  x  for  sudden 
temperature  change  in  T300/5208  laminate  having  a 
circular  cavity 


CONCLUDING  REMARKS 


The  major  contribution  of  this  work  may  not  lie  in  the  formulation  of  a  hy- 
grothermal  elasticity  theory  which  considers  coupling  between  mechanical  deforma¬ 
tion  and  diffusion  but  the  analytical  scheme  for  evaluating  the  physical  constants 
that  allowed  a  numerical  evaluation  of  the  results.  In  consistent  with  the  re- 
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cent  experimental  observations  [5],  the  applied  mechanical  stresses  can  strongly 
influence  the  moisture  distribution  within  composite  materials.  Such  observation 
obviously  suggests  a  strong  interaction  between  the  deformation  and  moisture  dif¬ 
fusion  process.  The  predictions  made  on  the  T300/5208  composite  laminate  with  a 
circular  cavity  show  that  stress  coupling  can  influence  the  results  both  quanti¬ 
tatively  and  qualitatively  by  a  wide  margin.  The  neglect  of  stress  coupling  leads 
to  unconservative  predictions  for  the  transient  moisture  and  temperature  boundary 
considered  in  this  investigation.  Such  findings,  however,  may  not  be  valid  in 
general  as  coupling  effects  are  sensitive  to  changes  in  the  combination  of  mois¬ 
ture,  temperature  and  stress  boundary  conditions.  The  influence  of  stress  coup¬ 
ling  is  ccuie.-6p&cyi^ic  and  should  be  investigated  by  analyses  and/or  experiments 
to  determine  the  extent  of  its  influence. 

The  present  work  should  be  regarded  only  as  a  preliminary  effort  towards  the 
development  of  a  fundamental  basis  for  understanding  the  continuum  parameters 
that  control  material  damage  in  the  presence  of  moisture  and  temperature  in  addi¬ 
tion  to  mechanical  stresses.  Moisture  enhanced  crack  growth  in  a  stressed  envi¬ 
ronment  is  one  of  the  least  understood  phenomena  in  fracture  mechanics.  There  is 
no  doubt  that  diffusion  and  mechanical  deformation  in  front  of  a  sharp  crack  under¬ 
go  strong  interaction.  Such  studies  are  presently  being  carried  out  in  the  In¬ 
stitute  of  Fracture  and  Solid  Mechanics  at  Lehigh  and  should  lead  to  a  better 
understanding  of  material  damage  enhanced  by  moisture  and  temperature. 


fl 
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